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Exciton effects on conjugated polymers and the resulting optical nonlinearities are investigated in soliton lattice 
states. We use the Su-Schrieffer-Heeger model with long-range Coulomb interactions for electrons. The third- 
harmonic generation (THG) at off-resonant frequencies is calculated as functions of the soliton concentration and 
the chain length of the polymer. The magnitude of the THG with 10 percent doping increases by the factor of 
about 10 2 from that of the neutral system. The large increase of the order two is common for several Coulomb 
interaction strengths, and is seen in open systems as well as in periodic systems. 
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The doping effects in conjugated polymers and their 
linear and nonlinear optical responses are fascinating re- 
search topics because of their importance in scientific 
interests as well as in developing technology. In the pre- 
vious paper, 1 ) we theoretically considered cxciton effects 
in the soliton lattice states of doped systems. There is 
one type of exciton in the undoped system with half- 
filled electronic states, 2 ' where the excited electron (hole) 
sits at the bottom of the conduction band (top of the 
valence band), which we called the intercontinuum ex- 
citon. 1 ' In the soliton lattice states of the doped Su- 
Schrieffer-Heeger (SSH) model for degenerate conjugated 
polymers, 3 -' there are small energy gaps between the soli- 
ton band and the continuum states, i.e., valence and 
conduction bands. The number of the types of excitons 
increases, and their presence effects structures of the op- 
tical spectra. A new exciton, which we have named the 
soliton-continuum exciton, 1 ' appears when the electron- 
hole excitation is considered between the soliton band 
and one of the continuum bands. 

We are conducting this study to investigate how 
large optical nonlinearities are obtained when conjugated 
polymers are doped with electrons or holes up to as much 
as 10 percent. The SSH model, 3 ' with the long range 
Coulomb interactions of the Ohno expression, 4 ' is solved 
with the Hartree-Fock (HF) approximation, and the ex- 
citation wavefunctions of electron-hole pairs are calcu- 
lated by the single-excitation configuration-interaction 
(single-CI) method. We calculated the off-resonant non- 
linear susceptibility as a guideline of the magnitude of the 
nonlincarity. We looked at the characters of the excita- 
tions at off-resonances resulting from the single-CI cal- 
culations. The third harmonic generation (THG) at the 
zero frequency, XthgW = X ( " 3 ' ) (3a->;a;,(jj,aj)| w= o, is con- 
sidered with changing of the chain length and the soliton 
concentration. We will demonstrate that the magnitude 
of the THG at the 10 percent doping increases by the 
factor of about 10 2 from that of the neutral system. 

The periodic systems simulate the reasonably long sys- 
tems where chain end effects are small. However, all of 
the polymers are not sufficiently long enough to neglect 
end effects, and shorter chains are present. The chain 
ends might have effects on the electronic structures, ex- 
citation properties and the resulting nonlinear optical 
responses of the polymers. In this paper, we discuss how 
the chain ends effect the large increase of optical non- 
linearity, compared with that of the periodic systems. 
We will demonstrate that the large increase of the THG 
upon doping is seen in open systems as well as in periodic 
systems. 

We use the SSH hamiltonian 3 ' with the Coulomb in- 
teractions, 



H = Hssh + Hi 
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The first term of eq. (1) is, 
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where t is the hopping integral of the system without the 
bond alternation; a is the electron-phonon coupling con- 
stant which changes the hopping integral linearly with 
respect to the bond variable y»; the bond variable is as- 
signed to the bond between the ith and (i + l)th sites 
(1 < i < N — 1, N being the system size); Ci j(T is an anni- 
hilation operator of the 7r-electron at the site i with spin 
a; the sum is taken over i within 1 < i < N — 1 due to 
the open boundary conditions; and the last term with the 
spring constant K is the harmonic energy of the classical 
spring simulating the tr-bond effects. The second term 
of eq. (1) is the long-range Coulomb interaction in the 
form of the Ohno potential: 4 ) 
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where n e \ is the number of 7r-electrons per site, r^j is the 
distance between the ith and jth sites, and 

W(r) = 1 (4) 

^(i/ur + (r/aVr 

is the Ohno potential. The quantity W(0) = U is the 
strength of the onsite interaction, V is the strength of 
the long range part, and a is the mean bond length. 

The model is treated by the HF approximation and the 
singlc-CI for the Coulomb potential. The bond variables 
are calculated by the adiabatic approximation. The self- 
consistent formalism was explained in our previous pa- 
per. 1 ' The electric field of light is parallel to the polymer 
chain which is along the x-axis. The optical absorption 
spectra are calculated by the formula 

J2E K P(uj-E K )(g\x\ K )( K \x\g), (5) 

where P(cu) = j/[7t(uj 2 +7 2 )] is the Lorentzian distribu- 
tion (7 is the width), E K is the electron- hole excitation 
energy, |k) is the Kth excitation, and |g) is the ground 
state. 

The THG is calculated with the conventional for- 
mula, 5 ~ 7 ) which is sometimes called the sum-over-states 
method. In order to demonstrate the magnitude of the 
THG, we use the number density of the CH unit, which 
is taken from trans-polyacetylene: 5.24 x 10 22 cm~ 3 . 8 ) We 
also use t = 1.8 eV in order to look at numerical data in 
the esu unit. We include a small imaginary part 77 in the 
denominator, which assumes a lifetime broadening and 
suppresses the height of the (^-function peaks. The THG 
at u) = does not sensitively depend on the choice of 77. 
This can be checked by varying the broadening. Here, 
we report the results with the value r\ — 0.02i. 

The system size is chosen as N — 80, 100, 120 when 
the electron number is even (it is varied from 7V i = 
N, N+2, N+A, N+6, N+8, N+10 to JV+12). We change 
Coulomb interaction parameters arbitrarily within a rea- 
sonable range in order to look at general properties of 
the optical nonlincarities of the soliton lattice systems. 
We take two combinations of the Coulomb parameters 
(U, V) — (2i, It) and (At, 2t) as representative cases. The 
other parameters, t = 1.8 eV, K = 21 eV/A 2 and a = 4.1 
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eV/A, are fixed in view of the general interests of this 
study. All the quantities of the energy dimensions are 
shown in units of t. 

Figure 1 shows the typical lattice configuration and 
excess-electron density distribution of open systems for 
N = 100, N eX = 102, U — At and V = 2t. Both quanti- 
ties are smoothed by removing small oscillations between 
even- and odd-number sites. There are two charged soli- 
tons due to the excess-electron number N c \ — N = 2. The 
two solitons are centered around the 20th and 80th lat- 
tice sites. They would be around the 25th and 75th sites 
in the periodic boundary system. The solitons are moved 
slightly closer to the chain ends. This is the end effect, 
due to the enhanced bond-alternation strengths near the 
1st and 100th sites. The enhanced bond variables pull 
the two solitons in the direction of the chain ends. These 
effects have already been seen in an earlier paper. 9 ) The 
excess-electron density shows that the doped charges ac- 
cumulate with the maxima at the soliton centers. 

Following, we calculate the optical spectra, linear ab- 
sorption and THG, and consider the exciton effects. Fig- 
ure 2 (a) shows the typical optical absorption spectrum 
at the 2 % soliton concentration for (U, V) = (At,2t). 
The broadening 7 = 0.05i is used. There are two main 
features around the energies 0.7t and I. At. The former 
originates from the soliton-continuum exciton, and the 
latter is from the intercontinuum exciton. The presence 
of the chain ends does not effect the excitation energies so 
much. However, the oscillator strengths of the soliton- 
continuum exciton become relatively larger than those 
in the periodic system. This is a new property which 
is found in the system with open boundaries. We have 
treated the periodic system with ring geometries, and 
the open system with geometries of linear chains. Thus, 
the expectation values of dipole operators are generally 
larger in the open system than in the periodic system. A 
similar fact has been discussed for the uniformly dimcr- 
ized system without solitons in rcf. 10. 

Figure 2 (b) displays the absolute value of the THG 
against the excitation energy u>. The abscissa is scaled by 
a factor of 3 so that the features in the THG are located 
at similar points to those in the abscissa of Fig. 2 (a). 
The large feature at approximately u> = Q.22t originates 
from the lowest excitation of the soliton-continuum ex- 
citon and the other feature at approximately u> = 0.26i 
originates from the higher excitations. The features from 
the intercontinuum exciton extend from u> = 0.48i to the 
higher energies. The soliton-continuum exciton gives rise 
to the large optical nonlincarity, as we looked in the lin- 
ear absorption. The oscillator strengths accumulate at 
the lowest excited state of the soliton-continuum exciton 
while the doping proceeds, as we have shown in ref. 1. 
This feature mainly contributes to the large optical non- 
linearity of the doped systems. The THG spectra, such 
as in Fig. 2 (b), are calculated for the three system sizes, 
N = 80, 100, 120, and for the soliton concentrations up 
to 10 %, in order to compare to the calculations for the 
system with periodic boundaries. 

Figures 3 (a) and 3 (b) display the variations of the ab- 
solute value of Xthg(°) for ( u > V ) = ( 2t > lt ) and ( 4 *> 2i )> 
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respectively. The plots are the numerical data: the open 
symbols are for the periodic systems, and the closed sym- 
bols are for the present calculations of the systems with 
chain ends. The dashed lines are guides for the plots of 
the periodic systems, showing the overall behavior for 
each system size. The finite system size effects 11 ) appear 
in these figures, ft is known that the THG is not size 
consistent, and spectral shapes depend on the system 
size when TV is as large as 100. n ) This fact is reflected in 
the separation of the plots of the three system sizes with 
periodic and open boundaries. 

The off-resonant THG near zero concentration in- 
creases very rapidly, but the THG still increases for a 
few percent up to a 10 % soliton concentration. The in- 
crease between the zero concentration and the 10 % con- 
centration is by the factor of approximately 100. This 
behavior is consistent for the two boundary conditions 
and the Coulomb interaction parameters. The main dif- 
ference between the two boundary conditions is that the 
THG with open boundaries is larger than that of the 
periodic system, when the concentration, system size, 
and the Coulomb parameters are the same. A similar 
property has been discussed in the calculations of the 
half- filled systems. 10 ) 

When we look at the enhancement factor in detail, it 
is discovered that the enhancement is particularly large 
for the open boundary systems with two solitons. This is 
seen for the three closed symbols near 2 % concentration 
in Figs. 3 (a) and 3 (b). These are the cases of the con- 
centrations of 2.5 %, 2.0 %, and 1.67 %, for N = 80, 100, 
and 120, respectively. Therefore, if the limit of the long 
chain is taken, there is a jump of the THG between the 
zero concentration and the infinitesimal concentration. 
This is different from the expected smooth behaviors of 
the periodic chains. The presence of the chain ends re- 
sults in the trapping of solitons near the ends, and thus 
the THG jumps suddenly at the zero concentration limit. 

The dramatic enhancement of the THG in the case of 
two solitons is clear when we consider the ratio of the 
THG value of the open system with respect to that of 
the periodic system. The calculations have been done 
for (U,V) = (At, 2t), and the ratio is shown in Fig. 4. 
The numerical data are shown by the triangles (N = 80), 
circles (N — 100), and squares (N = 120). Except for 
the three plots near the concentration of 2 %, the ratio 
is constant at approximately 3. It becomes greater than 
10 for the two soliton systems. Thus, we have found 
that the case with two solitons is special, mainly due to 
the chain end effects. It should be noted that the large 
THG due to the presence of chain ends could be used 
as a tool for increasing nonlinear optical responses ex- 
perimentally. The reason for the large enhancement can 
be understood as follows. In two soliton solutions with 
the open boundary condition, the intersoliton distance 
becomes longer than in the system with periodic bound- 
aries. This gives rise to the larger expectation values of 
the dipole moment operators, and thus we obtain the 
greatly enhanced THG in the system with two solitons. 

In summary, we have considered the off-resonant non- 
linear susceptibility as a guideline of the strength of the 
nonlincarity in the doped conjugated polymers. We have 
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calculated the off-resonant THG with various system 
sizes and the soliton concentration for the chains with 
open boundaries. We have shown that the magnitude of 
the THG at 10 percent doping increases by the factor of 
10 2 from that of the neutral system. The large increase 
of the order two is common for several Coulomb interac- 
tion strengths, and is seen in the open systems as well as 
in the periodic systems. 
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Figure captions 

Fig. 1. (a) The bond alternation order parameter, (—l) n (j/ n +i—j/ n )/2, 
and (b) the excess-electron density, (p n -l + 2pn + pn+l)/4. The 
parameters are U = it, V = 2t, N = 100, and TV e] = 102. See 
the text for the other parameters. 

Fig. 2. (a) The optical absorption spectrum and (b) the absolute 
value of the THG, for the open system with the size TV = 100, 
the electron number TV cl = 102, and (U, V) = (4t,2t). The 
broadening 7 = 0.05t is used in (a), and r\ = 0.02t is used in (b). 
The absorption is shown in arbitrary units, and the nonlinear 
optical response is in esu units. 

Fig. 3. The absolute value of the THG at ui = v.s. the soli- 
ton concentration for (a) (U, V) = (2t, It) and (b) (4t,2t). The 
numerical data are shown by the triangles (N = 80), circles 
(N = 100), and squares (TV = 120). The data of the system with 
the periodic boundary condition are shown by the open symbols. 
And, the data with the open boundaries are shown by the closed 
symbols. The dashed lines are guides. 

Fig. 4. The ratio of the absolute values of the THG between 
the open and periodic boundary conditions, shown against 
the soliton concentration. The Coulomb parameters arc 
(U, V) = (At, 2t). The numerical data are shown by the triangles 
(TV = 80), circles (TV = 100), and squares (TV = 120). 
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